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Abstract. The effects of tetragonal strain on electronic and magnetic properties of 
strontimn-doped lanthanum manganite, La2/3Sr]^/3Mn03 (LSMO), are investigated by 
means of density-functional methods. As far as the structural properties are concerned, 
the comparison between theory and experiments for LSMO strained on the most 
commonly used substrates, shows an overall good agreement: the slight overestimate 
(at most of 1-1.5 %) for the equilibrium out-of-plane lattice constants points to possible 
defects in real samples. The inclusion of a Hubbard-like contribution on the Mn d 
states, according to the so-called "LSDA+U" approach, is rather ineffective from the 
structural point of view, but much more important from the electronic and magnetic 
point of view. In particular, full half-metallicity, which is missed within a bare density- 
functional approach, is recovered within LSDA-I-U, in agreement with experiments. 
Moreover, the half-metallic behavior, particularly relevant for spin-injection purposes, 
is independent on the chosen substrate and is achieved for all the considered in-plane 
lattice constants. More generally, strain effects are not seen to crucially affect the 
electronic structure: within the considered tetragonalization range, the minority gap 
is only slightly {i.e. by about 0.1-0.2 eV) affected by a tensile or compressive strain. 
Nevertheless, we show that the growth on a smaller in-plane lattice constant can 
stabilize the out-of-plane vs in-plane eg orbital and significatively change their relative 
occupancy. Since eg orbitals are key quantities for the double-exchange mechanism, 
strain effects are confirmed to be crucial for the resulting magnetic coupling. 
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1. Introduction 

Strontium-doped lanthanum manganites, Lai-xSr^^MnOs, have been heavily studied 
in the past few decades [H 12 E], due to their fascinating properties and intriguing 
apphcations. Nevertheless, most of the properties are composition-dependent and the 
most exotic ones, such as colossal magneto-resistance and/or half-metallicity, appear 
only when 0.17< x <0.5[2]. In the case of the optimized composition (x ~l/3), 
Lao.eeSro.ssMnOs (denoted as LSMO hereafter) shows a spin polarization >95%|4| and 
a high Curie ferromagnetic temperature {Tc ~ 370 K); therefore, this compound appear 
as a promising candidate in the spin injection framework, especially since it is envisaged 
that LSMO-based devices could operate at room temperature. 

The interactions between electrons and lattice distortions are well known to play a 
significant role in the physics of the compound: Mn Cg electrons, which are responsible 
for the rich variety of attractive properties shown by manganites, are coupled to the 
lattice through Jahn- Teller effects [Sj. In particular, MnOe octahedral tetragonal 
distortions can be induced by the lattice mismatch between the film and the substrate, 
to which magnetic anisotropy and magnetoresistance effects of LSMO films were found 
to be extremely sensitive [H]. In particular, the magnetization easy axis was shown to 
move from in-plane to out-of-plane upon moving from tensile to compressive strain 
[3 IE] . Moreover, it is mandatory for industrial applications to achieve a good control 
over the film growth, which in turn requires a careful choice of appropriate substrates 
to get optimized LSMO films. High-performance epitaxial manganites thin films 
have been successfully grown on [001]-oriented oxides, such as LaAlOa fLAO)P ITU], 
Lao.3Sro.7Alo.35Tao.35O9 (LSAT)^6j, NdGa03 (NGO)^, and SrTiOg (ST0)|11 HH Cg. 
Strained manganites epilayers often show a different behavior with respect to the 
bulk. In particular, for LSMO, tensile strain was shown to reduce the Curie ordering 
temperature, Tc, and this was successfully explained in terms of double exchange: 
the increase of the in-plane Mn-0 bond length leads to a reduction of the hopping 
term between Mn^^-Mn^^, therefore reducing Tc (see below) |5]. However, for other 
manganites, such as Lai_a;Baa;Mn03, a tensile strain was found to enhance Tc [HI- 
Therefore, a thorough understanding has not yet been achieved, due to contradictory 
results. Possible reasons might include the still poor control over the samples growth and 
stoichiometry. Moreover, one should keep in mind that the delicate interplay between 
charge, spin and orbital ordering can be very easily altered in manganites and the 
subtleties of the relevant physical effects (such as magnetoresistance and anisotropy) 
can lead to strikingly different results induced by generally tiny changes. 

Despite the large abundance of experimental studies devoted to substrate-induced 
strain effects on the LSMO magnetic properties [HI GDI and ab-initio theoretical works 
focused on LSMO,jini CH CHI HEl IHI a deep investigation of strain effects on 
the electronic and magnetic properties from first-principles (both with and without 
additional correlation effect beyond standard exchange-correlation functionals) is lacking 
and will therefore be the focus of the present work. We choose four typical substrates 
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whose mismatches with LSMO films induce different types of lattice strain: LAO, for 
which the film undergoes a relatively large in-plane compressive strains; LSAT and 
NGO, for which the film undergoes a very small compressive strain; STO, for which the 
film experiences a small tensile strain. The paper is organized as follows: in Sec|21we 
report some technicalities (in terms of structural and computational details); in SecOlwe 
focus on the structural equilibrium properties of LSMO strained on different substrates. 
Then, we move in Sec|3]to the discussion of the effects of strain and correlation on the 
LSMO electronic and magnetic properties, in terms of density of states and magnetic 
moments. Finally, we summarize our conclusions in SeclHl 

2. Crystal structure and computational details 

Our calculations are performed following the generalized gradient approximation 
(GGA)PH] to the exchange-correlation potential within the framework of density 
functional theory (DFT), using the all-electron full-potential linearized augmented 
plane- wave (FLAPW)^n] method as implemented in the FLEUR-code^|. In order 
to take into account correlation effects beyond the local spin density approximation 
(LSDA) or the GGA, we used the so called "LSDA+U" (or "GGA+U", in this case) 
scheme pTj in the "atomic limit" approximation, as implemented in FLEUR[221; this 
formalism includes a Hubbard-like potential acting on the Mn 3d states that enhances 
the tendency towards localization which erroneously lacks in bare GGA. 

The unit cell is divided into non-overlapping spheres and interstitial region. The 
muffin-tin (MT) radii are set to 2.5 a.u. for Sr and La, 2.0 a.u. for Mn and 1.5 a.u. for 
O, respectively. The wave function expansion in the interstitial region was carried up to 
kmax = 3.8 a.u.~^ which limits the number of basis functions to be 600 or so. We treated 
La 5s, 5p; Sr 4s, 4p as local-orbital states^23j, whereas La 4/ were treated as valence 
states. The energy position of these latter states (resulting to be at ~2 eV above the 
Fermi level, see below) can be questionable and might be an artifact of DFT-FLAPW; 
however, tests made applying a large Coulomb parameter (U=10 eV) on La 4/ states 
show that the effect on the physical quantities we are mainly interested in [i.e. total 
energies and density of occupied states) is neghgible. 

In order to sample the irreducible Brillouin zone, 72 k-points were used, according 
to the 12x12x4 Monkhorst-Pack shell|21]; the convergency with respect to this parameter 
was accurately tested. Within the GGA+U formalism, we considered different values of 
U (such as 2 and 3 eV), which were found to give the better agreement of DFT results 
with recent photoemission measurements. [221 The J value was set to J = 0.7 eV. 

In order to simulate LSMO, we used a 15-atoms unit cell which can be schematically 
represented as a [001]-ordered (SrMnOs)!/ (LaMn03)2 superlattice (see FigH}, similar to 
previous works^S]- With this choice of the unit cell, there are two inequivalent Mn: Mui 
has a LaO plane on one side and a SrO plane on the other side, whereas Mn2 has LaO 
layers on both parts. Possible MnOg octahedral tiltings were neglected. Effects induced 
by this approximation are quantitatively hard to predict; moreover, experiments suggest 
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Figure 1. Schematic unit cell (Bravais lattice shown by thin solid lines) for LSMO: 
grey, blue, red and black spheres denote Sr, O, Mn and La atoms, respectively. The 
vertical axis shows the [001] direction. 0-0 bonds are shown via blue rods. 

the tilting angle not to be strongly affected by different strain conditions. The in- 
plane lattice constants a are chosen according to the LAO, NGO, LSAT and STO bulk 
experimental values^6j this is of course an approximation, and the theoretical substrate 
lattice constants might be used, as well. In some cases {i.e. when the disagreement 
between the theoretical and experimental lattice constants for the substrate or for the 
epilayer is rather high), one or the other choice can lead to different strain signs and 
qualitatively different behavior. However, in the present case, when considering LSMO 
grown on different substrates, the choice of the experimental substrate lattice constants 
does not change the strain sign with respect to the experimentally realized situation. 
In fact, the LSMO equilibrium calculated (see below) and experimental "cubic" lattice 
constant are a^J^ = 3.90 Ang and Cg^^* = 3.874 Ang, respectively. Therefore, in both 
simulations and experiments, the epilayer is under tensile strain for STO and under 
compressive strain in the other cases. The out-of-plane lattice parameters c have been 
set to i) experimental values for LSMO thin films jHI grown on different substrates and 
a) minimized according to the calculated total energies (see below). To test the effect of 
structural optimization, we relaxed the atomic positions of LSMO strained on STO until 
the forces were less than 0.002 Hartree/(a.u.). We verified that the electronic structure is 
negligibly affected by structural optimization of internal atomic coordinates; therefore, 
we will focus in the following on the unrelaxed structures. In all the simulations, 
the three Mn atoms in the unit cell were aligned ferromagnetically, accordingi2j to 
experiments P7] . 




Figure 2. a) Out-of-plane lattice constants (denoted as c) vs in-plane lattice constant 
(denoted as a): experimental values, GGA and GGA+U results are shown as blue 
stars, black filled circles and red filled diamonds, respectively. Linear fits to each set 
of data are also shown. The total energy curves as a function of c are shown (with 
black symbols and solid lines) for different substrates in panels b), c), d) and e) for 
LSMO LAO , LSMO jvGO , LSMOlsat and LSMOsto , respectively. In b) and e), the 
red symbols and dashed line show the results of GGA+U calculations (U=2 eV). Axes 
labels in panels b), c) and e) (not reported for clarity in plotting) are the same as d). 



3. Structural equilibrium properties 

In FigI21we show the total-energy vs out-of-plane lattice constant for different substrates 
(see panels b), c), d) and e) for LSMOlao , LSMOatgo , LSMOL^^iT and LSMO^ro j 
respectively); we performed a parabolic fit of the total energies and obtained the 
equilibrium c as the parabola minimum for each different substrate. The resulting c 
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values are then plotted in FigI21 function of the substrate lattice constants and 

compared with corresponding experimental values. Similarly, in Tabled we report the 
in-plane lattice constants (a) for the different chosen substrates, as reported in Ref. [Oj, 
along with the corresponding experimental and theoretical values for the out-of-plane 
lattice constant (c), the strain components (en), the c/a ratio as well as the resulting 
volume {V). The LSMO equilibrium [i.e. unstrained) lattice constant was estimated as 
the weighted average between the experimental lattice parameter of the constituents, 
LaMnOs and SrMnOs and is evaluated as aeq = 3.874 A.[TT] 

First of all, as expected from elasticity theory, the c parameter decreases as 
the strain becomes more and more tensile. However, both from the theoretical and 
experimental point of view, the lack of a constant volume shows that the material 
does not behave as an ideal elastic medium: the more compressive (tensile) is the 
strain, the smaller (larger) is the volume. Moreover, we see from FigEl a) that in the 
considered range of substrate lattice constants, the c vs a relation is well described using 
a linear fit, both for theoretical as well as experimental values. However, the slope of 
the experimental results is larger than the theoretical one. In particular, GGA slightly 
overestimates the c parameter in all cases (at most by 1.3%), except for LSMO^ao- This 
can be interpreted as follows: on one hand, the slight disagreement between experimental 
and theoretical c values might be due to non-perfect experimental samples. In fact, it 
is well known that twin-like defects or oxygen lack - or excess - can often occur 
as a result of the different growth conditions. The experimental films might therefore 
show a non-perfect stoichoimetry or a non-high-quality crystallinity (at variance with 
our simulated perfect crystal), therefore resulting in slightly different lattice constants. 
On the other hand, we have to keep in mind that GGA is known to describe rather well 
[i.e. within 1-2 %) the structural properties of magnetic systems (generally better than 
LSDA), though showing, in some cases a slight overestimation of the lattice constants. 
This could be the case, although the disagreement with experiments is within the 
generally acknowledged error range. Therefore, in order to check whether the error 
is due to the neglect of correlation effects, we repeated the calculations for the LAO and 
STO substrates (see red diamonds and dashed lines in FigOb) and e), respectively), 
using a GGA-I-U approach with U=2 eV. In both cases, we found a slight increase of the 
c parameter, which brings LAO in excellent agreement with experiments, but STO to 
an even worse value. Our results are perfectly consistent with previous studies on other 
magnetic oxides j2H]: the addition of a Hubbard-like potential to bare DFT generally 
leads to only slight changes in the structural propertes (and in particular to a larger 
volume, consistently with our results), though it is generally much more "effective" 
from the electronic structure point of view, bringing about a much closer agreement 
with spectroscopic properties, compared to DFT (see below). 
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Table 1. Experimental and theoretical structural parameters for LSMO films 
grown on different substrates: in-plane lattice constants (a, in A), out-of-plane lattice 
constants (c, in A), c/a ratio and unit-cell volume {V, in A'^). Experimental values 
are taken from Ref. 6 . In parenthesis, we report the in-plane and out-of-plane strain 
components defined as e^x = {0' — a-eq)/o.eq and Ezz = ic— a^q) j a^q^ respectively. Strains 
are evaluated with respect to the ideally unstrained experimental cubic lattice constant 
(Oe, = Ceg = 3.874 A). 



Substrate 


a {(-xx) 


Expt. Th. 


c/a 
Expt. Th. 


V 

Expt. Th. 


LAO 
NGO 
LSAT 
STO 


3.793 (-2.1%) 
3.861 (-0.3%) 
3.868 (-0.1%) 
3.905 (0.8%) 


3.991 (3.0%) 3.971 
3.902 (0.7%) 3.927 
3.892 (0.5%) 3.924 
3.850 (-0.6%) 3.899 


1.052 1.048 
1.011 1.017 
1.006 1.014 
0.986 0.999 


57.42 57.19 
58.17 58.54 
58.23 58.71 
58.71 59.47 



4. Electronic and magnetic structure 

Before starting the discussion of strain-induced effects on tlie magnetic and electronic 
properties, we briefly comment on one of tlie systems of interests, namely LSMO strained 
on LAO, in terms of density of states and magnetic moments. Next, we will focus on 
the comparison between this same system and others under different strain conditions 
{i.e. different substrate lattice constants), according to the experimental equilibrium 
structure; finally, we will focus on the systems on a fixed in-plane lattice constant and 
study the effects of the c elongation. 

4.1. Electronic and magnetic structure of LSMO strained on LAO: hybridization and 
correlation effects 

We show in Fig. El the LSMO^yio density of states projected (PDOS) on the relevant 
atoms: Mn and O. As well known. La and Sr only donate their electrons to the electronic 
system, the bonding being dictated essentially by the Mn-0 interaction. In fact, as 
clearly visible in FiglHl there is a strong hybridization between Mn d and O p states, 
especially in the majority occupied states. The peak at binding energies at -1.2 eV and 
-1.8 eV show the t2g orbitals (not so strongly mixed), whereas the remaining part of 
the Mn PDOS shows the Cg contributions, much broader due to hybridization. As can 
be also seen from the PDOS, and will be better clarified in the following, the system is 
"nearly half- metallic" , i.e. it shows a gap of ~ 1-2 eV (see below), but the conduction 
band minimum (CBM) is just below the Fermi level {Ep). Therefore, the system is 
not strictly 100% spin-polarized at -E^.j^n] The situation changes upon introduction 
of the Coulomb correlations: as shown in FigOl in the case of U = 3 eV, the system 
becomes fully half-metallic, consistently with experiments suggesting a degree of spin- 
polarization larger than 95%|lj. As expected, the occupied (unoccupied) d states are 
shifted towards higher binding (higher excited) energies by about 1 eV {i.e. of the order 
of (U-J)/2) and, due to the strong hybridization, O p states are also affected (see FigOl 




-8 -6 -4 -2 2 -8 -6 -4 -2 2 -8 -6 -4 -2 2 
E (eV) E (eV) E (eV) 

Figure 3. PDOS for a) Mni, b) Mna, c) Oi, d) O2, e) O3 and f) O4 for LSMOlao ■ 
Black solid and red dashed lines show the results of GGA and GGA+U calculations, 
respectively (U= 3 eV). Atomic labels are consistent with Fig^ 

c), d), e) and f)). Therefore, since within bare GGA half-metalhcity is broken because 
the tails of minority Mn d states cross Ep, the Hubbard-induced shift towards higher 
energies results in a fully half-metallic character for LSMOlao ■ 

As for the magnetic moments, we recall [201 ^| that the two Mn atoms are 
surrounded by different cationic environments: Mni has a LaO plane on one side and a 
SrO plane on the other side, whereas Mn2 has LaO layers on both parts. As a result, the 
magnetic moment of Mn2 (3.26 fis) is closer to pure LaMnOa in similar strain conditions 
{i.e. of the order of 3.3 /ib), whereas the moment of Mui (3.11 fis) is, as expected, lower 
than Mn2, but somewhat larger than the average of LaMnOa and SrMnOs (this latter 
showing a muffin-tin moment of about 2.7 fis)- The changes in the electronic structure 
induced by the presence of additional correlations (and in particular the shift of Mn d 
states towards higher binding energies) affect in turn the magnetic moments: whereas, 
in the bare GGA case, the total magnetic moment is 10.89 /i^, this increases up to 
11.0 /is, when the GGA+U approach is used. The integer total magnetic moment 
is consistent with the fully half-metallic character previously outlined. Similarly, the 
atomic moments increase up to 3.29 fiB and 3.49 fiB for Mui and Mn2, respectively. 
However, neither within GGA nor within GGA+U, the picture of charge ordering and 
mixed Mn^"*" and Mn^+ valencies is valid within the present DFT-based formalism and 
with the present considered unit cell, as already outlined in previous papers^Hj- Finally, 
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Figure 4. TDOS for LSMO grown on LAO (black solid line), LSAT (blue dashed line) 
and STO (red dot-dashed line)with GGA method. The zero of the energy scale is set 
to the Fermi level. Majority (minority) spin contributions are shown in the positive 
(negative) y-axis. For clarity, the TDOS in the case of LSMOsto and LSMOlsat are 
arbitrarily shifted along the y-axis. Vertical dotted lines are guide to the eyes. 

we remark that the magnetic moment is about 3.66 yU^/Mn, in rather good agreement 
with experimentally reported values of 3.7 /i^ [TIJ. 

Strain effects on the electronic and magnetic structure 

Let us now focus on the structures grown on different substrates at the experimental 
a and c lattice constants (see Table Q). In Fig^we report the total density of states 
(TDOS) of the LSMO under tensile and compressive strains. The TDOS for LSMO 
grown on NGO looks very similar to LSMO grown on LSAT and is therefore not shown. 

The comparison shows that the TDOS all look rather similar. However, one 
can note that, in going from compressive to tensile strain, the band structure (both 
in the majority and minority spin channels) shifts towards lower binding energies 
(by ~few tenths of eV for the considered strains of few percents). Furthermore, in 
none of the considered substrates, the system shows full half-metallicity; rather, all 
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of them are "nearly" half-metallic. In order to further investigate the role of strain 
on the half-metallic behavior, we report in Fig. |31 the energy position (taking Ep as 
reference) of the valence band maximum (VBM) and of the CBM in the minority spin 
channel as a function of the substrate lattice constant. The plots show that in all the 
considered systems, the CBM is always lower than Ep, therefore hindering the full half- 
metallicity. However, there is a slight tendency towards half-metallicity upon applying 
a tensile strain: LSMO grown on STO is closer to half-metallicity than LSMO grown on 
LAO. Furthermore, the minority band-gap slightly increases with the substrate lattice 
constant. 

Let us now focus on the effects of correlation and of strain on the band-gap and of 
the relative position of Ep with respect to the minority band edges. Upon introduction 
of U and already for a small value U = 2 eV, the Fermi level enters the gap. As expected, 
the gap (here defined as the difference between the CBM and the VBM, irrespective 
of the position of Ep) increases with applying progressively increasing U parameters: 
within the range of considered strains [i.e. in going from LSMOpAo to LSMOsto )? it 
is of the order of 1.2-1.4 eV within bare GGA, of 1.6-1.8 eV within GGA+U (U = 2 
eV) and of 1.9-2.1 eV within GGA+U (U = 3 eV). On the other hand, no significative 
changes are observed upon applying strains (at least in the experimentally accessible 
range): the strain-induced effects on the gap-magnitude amount at most to 0.1-0.2 eV. 

It is of further interest to investigate the effects on the strain of the magnetic 
properties, in terms of magnetic moments. In going from LSMO^ao to LSMO^^o ? the 
Mni and Mn2 moments range from 3.11 /is to 3.15 fj,B and from 3.26 /Xb to 3.29 /is, 
respectively. In parallel, the total moment ranges from 10.89 hb to 10.94 /z^. Therefore, 
consistently with the slight changes in the electronic structure, the magnetic moments 
only slightly increase upon increasing the substrate lattice constant. 

Since the physics of LSMO crucially depends on the Cg orbitals, whose occupancy 
and itineracy strongly affect the resulting ferromagnetism through double exchange, 
we further investigate the Mn d PDOS, further resolved by m character. We recall 
that in a cubic symmetry, the Cg level is twofold degenerate; however, upon tetragonal 
strain, the degeneracy is lifted and the — y'^ in-plane orbital becomes different from 
the 3z'^ — out-of-plane orbital. Therefore two competing mechanisms are induced 
by strain: i) on one hand, according to double exchange, Tc is proportional to the 
transfer integral t^, describing the Mn^^-Mn^"*" hopping through intervening oxygen 
atoms. In turn, to oc (ig ^'^ where do is the in-plane Mn-0 bond length. Therefore, a 
tensile strain, which increases the in-plane do, should reduce to, and, as a consequence, 
Tc', ii) on the other hand, the relative occupancy of the — and 3z^ — orbitals 
also depend on strain: in particular, for a tensile strain, the occupancy of the in-plane 
orbital should be higher than that of the out-of-plane orbital. In parallel, the — y"^ 
has a larger transfer intensity: the stabilization and related increased occupation of this 
orbital therefore enhances electron hopping and, therefore, Tq- The investigation of the 
strain-dependence of ferromagnetic coupling and of Tc goes beyond the scope of the 
present work. However, we provide here some hints on the occupation of the — y^ 
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Figure 5. CBM (red empty circles) and VBM (black filled circles) positions in the 
minority spin channel with respect to Ep (set to zero) as a function of the in-plane 
lattice constant (labels on top show the experimentally considered substrates) . Panels 
a), b) and c) show the results of bare GGA, GGA+U (U=2 eV) and GGA+U (U=3 
eV) approaches, respectively. 
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and 3z^ — r^ orbitals, which could serve as a starting point for further investigations on 
the magnetism in LSMO under strain conditions. 

In Figini we show the orbital-resolved DOS projected on the two Cg orbitals 
for LSMO5TO and LSMO LAO (both with and without the inclusion of the Hubbard 
parameter U) for the two inequivalent Mni and Mn2 atoms. We also show the 
integral of the PDOS, defined, for example for the in-plane Cg orbital of Mui, as: 
= f PDOS^2_y2{E') dE'. This is particular meaningful for E = Ep, where 
the integral coincides with the occupation of the orbital. Note that the t2g orbitals as well 
are affected by strain and the three-fold degeneracy lifted. However, their occupation 
is not modified upon tetragonal distortion: in fact, their PDOS (not reported) shows 
a peak at in the [-1.5:-l]-eV range and becomes negligible above ~-0.5 eV, resulting in 
full occupation of the majority subband contributions. On the other hand, minority t2g 
states are completely unoccupied (as already mentioned, their tails slightly cross Ep but 
this is likely to be an artifact of DFT and fully zero-occupation, i. e. half-metallicity, can 
be restored by applying an additional Coulomb correlation). It is evident from FigEla) 
that, even in LSMOsro; showing an in-plane tensile strain of +0.8%, the degeneracy is 
lifted, but the consequences are still small: as expected, the in-plane Cg is more occupied, 
but only by about 0.01 electrons more than in the out-of-plane Cg case. This picture is 
kept more or less unaltered even upon application of U (see Fig. IHlb)). On the other 
hand, if we look at LSMOlao , where the strain is markedly compressive, there is a 
considerable difference between the occupations of the Cg orbitals and the situation is 
reversed with respect to tensile strain: the Sz"^ — is about 15% more occupied than 
— y"^ for Mui within bare GGA (cfr FiglHlc)) and this same difference increases to 
20% in the case of GGA+U (cfr FiglBl d)) . Similarly and even more markedly, the 
occupation of the out-of-plane Cg orbital is 20% higher for Mn2 and further increases 
within GGA+U (cfr FiglHl e) and f)). This is expected to have consequences on the 
ferromagnetism (although it is hard to give a quantitative estimate on the reduction of 
Tq induced by the different occupations): in particular, upon compressive strain, the 
resulting lower transfer integral for the more-occupied — orbital should suppress 
the in-plane carrier density and, as a consequence, double-exchange hopping and Tc, 
whereas, as already pointed out, a smaller in-plane MnO bond- length should lead to an 
opposite tendency. 

Finally, we briefly comment on the effects of the Jahn- Teller Qs-like distortions |31j 
{i.e. octahedral elongation) for LSMOlao , LSMOatgo and LSMOl^io , in particular 
on the magnetic moments. In FigjTI we plot the magnetic moments (on the two Mn 
atoms and total moment per unit cell in panels b), c) and a), respectively) as a function 
of the c lattice constant for LSMO on the LAO, NGO and STO substrates, ; we also 
show the experimental equilibrium structures using filled symbols. As previously noted, 
the situation in LSMO lsat is very similar to LSMOatgo and is therefore not shown. It 
is worthwhile remarking that, whereas moments in LSMO^^o are basically unaltered 
by the out-of-plane lattice constants, LSMOatgo shows slightly larger strain- induced 
effects and in LSMO^^io the dependence on c is rather significative. This is consistent 




E (eV) E (eV) E (eV) 

Figure 6. Mni PDOS of — y"^ (black bold solid) and 3z^ — (red bold dashed) for 
LSMO on a) STO within GGA, b) STO within GGA+U (U=3 eV), c) LAO within 
GGA and d) LAO within GGA+U. The Mn2 PDOS of - (black bold solid) and 
3z2 - r2 (red bold dashed) for LSMO on e) LAO within GGA and f) LAO within 
GGA+U (U=3 eV). Also shown are the integrals: /*|^ 2(-B) (black thin solid) and 
/3f2"_^2(£^) (red thin dashed). 

with the related electronic structure: whereas LSMO on STO is always very close to 
half-met allicity for all the considered c values, this condition is by far "less satisfied" 
in LSMO grown on LAO. This has a relatively big effect on the magnetic moments: 
especially for lower c values, there is an appreciable density of states coming from the 
minority conduction bands that cross Ep, therefore becoming occupied and lowering 
the total as well as Mn-projected magnetic moment. 

5. Conclusions 

First-principles calculations have been performed on LSMO grown on different 
experimentally common substrates, in order to highlight the effects of uniaxial or biaxial 
strain on the electronic and magnetic structure. Our results show that: i) GGA 
does not exactly reproduce the experimental structural parameters, such as the out- 
of-plane lattice constant, and the inclusion of correlation effects according the GGA+U 
formalism does not significatively alter the situation. However, the discrepancy (at most 
of the order of 1.3 %) falls within the generally acknowledged error range of density 
functional calculations and might also be due to defects in the experimental samples; ii) 
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Figure 7. a) Total, b) Mni and c) Mn2 magnetic moments (in Bohr magnetons) as 
a function of the out-of-plane lattice constant (in A) for LSMOlao (black circles), 
LSMO jvGO (red diamond) and LSMOsto (blue squares). Filled symbols show the 
values obtained for the experimental equilibrium structure. 



in none of the considered substrates, we obtain full half-metallicity within bare GGA. 
However, even the inclusion of a small Hubbard parameter U ~ 2-3 eV - consistent 
with recent experiments - turns the system in a 100% spin-polarized density of states at 
the Fermi level, in agreement with experiments and for all the considered substrates . 
Moreover, within bare DFT, LSMOsto is closer to half-metallicity than LSMOlao '■ the 
minority t2g states, whose tails cross the Fermi level preventing full half-metallicity, have 
increasingly higher excited-state energies as the tensile strain grows. In summary, strain 
effects (at least in the experimentally accessible range) on the relevant electronic and 
magnetic properties, such as the half-metallic gap and the magnetic moments, do not 
significatively change the overall picture. However, subtle effects such as the occupation 
of the in-plane vs out-of-plane Cg orbital, that are quantitatively estimated in the case 
of LSMOlao and LSMOsto ? can change dramatically the magnetic couphng. These 
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results call for further investigations on the strain-induced effects on the ferromagnetic 
exchange constants. 
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